Spin-dipole induced lifetime of the least-bound 5 E+ state of He(2 3 Si)+He(2 3 Si) 
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The properties of the least-bound vibrational level (y = 14) of the E+ state formed during 
the ultracold collision of two spin-polarized metastable 2 3 Si helium atoms are crucial to studies of 
photoassociation spectroscopy of metastable helium. We report a calculation of the autoionization 
lifetime t 9 of this state induced by spin-dipole coupling of the 5 E+ state to the 1 Ti^ state from which 
Penning and associative ionization processes are highly probable. We find r 9 ss 150 /is, significantly 
larger than the recent experimental estimates of (4 — 5) (is. 
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I. INTRODUCTION 



Two-photon photoassociation spectroscopy of ultra- 
cold spin-polarized metastable 2 3 Si helium atoms has 
recently been used [l| to accurately measure the bind- 
ing energy of the least bound vibrational level (v = 14) 
of the 5 E+ state formed during the collision of two atoms. 
The measured binding energy Ef,(v — 14) = 91.35 ± 0.06 
MHz, combined with the new ab initio 5 £+ potential of 
Przybytek and Jeziorski [2|, yielded the high precision 
value a = 7.512 ± 0.005 nm for the s-wave scattering 
length. 

Knowledge of the decay width 7 9 of this bound state 
is important to the analysis of several experimental tech- 
niques used to study photoassociation in metastable he- 
lium. In two-photon photoassociation 0, yj, for the 
case where one laser is detuned far from resonance, the 
width of the narrow stimulated Raman line is deter- 
mined largely by ~f g , and, in dark resonance experi- 
ments, the incomplete suppression of the photoassocia- 
tion signal at the resonance minimum is due to 7 g . From 
these two experiments Moal et al 1] obtain the estimate 
0.05 < 7 9 /27r < 0.3 MHz, corresponding to a lifetime T g 
between 0.5 and 3 /us. The decay width is also needed 
in the interpretation of one-laser photoassociation exper- 
iments using calorimetric detection [3, la, 13 as decay of 
the v = 14 bound state contributes to the heating of the 
trapped atoms. 

Koelemeij and Leduc [3j have given a theoretical es- 
timate of T g based upon the semiclassical vibration fre- 
quency of the v = 14 level and inelastic decay rates due to 
spin relaxation and relaxation-induced ionization 0, Q ■ 
They obtain r g = (4 — 5) /is, in good agreement with 
the measured linewidths. However, these inelastic decay 
rates were calculated for initial scattering states in the 
5 E+ potential, not for a bound vibrational state. In this 
paper we shall calculate the decay of the v = 14 bound 
state due to relaxation-induced ionization and argue that 
there is no contribution from spin-relaxation. 



II. PERTURBATION THEORY 

A collision between two atoms with spin quantum num- 
bers Si and 6*2 involves the total spin state \(SiS2)SMs) 
where Ms is the projection of S onto the molecular axis. 
For metastable helium atoms Si = S2 = 1 and Penning 
and associative ionization processes are forbidden from 
the spin-polarized | (11)22) quintet state. However, the 
spin-dipole interaction between the metastable helium 
atoms couples the spin-polarized \SMg) = 1 22} state to 
the 1 21), 1 20) and |00) electronic spin states. As Pen- 
ning and associative ionization from the S = singlet 
state are highly probable, the spin-dipole induced cou- 
pling results in a finite lifetime of the 1 22} state. This 
spin-dipole interaction can be treated as a perturbation 
and the method presented here is a adaptation of that 
given in [9j. 

The total Hamiltonian for the relative motion is 

H = H Q + H sd , (1) 

where the Hamiltonian Hq for the unperturbed system is 
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2M 



V AS 



(r). 



(2) 



Here M is the reduced mass of the system and the adia- 
batic potential for the molecular state 2S + 1 A formed dur- 
ing the collision is denoted by V\s {f) where A is the pro- 
jection of the total electronic orbital angular momentum 
along the molecular axis. The spin-dipole interaction is 





H k a = 



h r 3 



[3{S 1 -r)(S 2 -f)-S l -S 2 ] 



(3) 



where Si and S2 are the electronic-spin operators for the 
two atoms, r = r/r is a unit vector directed along the 
internuclear axis and 

■*(*L) 2 E h al, 

where a is the fine structure constant, 

radius, (/i e //xs) is the electron magnetic moment to Bohr 

magneton ratio and Eh = a 2 m e c 2 is the Hartree energy. 
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(4) 
ao is the Bohr 



The eigenstates \a) of Ho can be written 
\a)=r- 1 R va {r) |* ) 



(5) 



where a = (v,a), and a = {T, S, Ms, l,m} labels the 
channel states 



\* a ) = \rSMs)\lm) 



(6) 



involving the relative motion states \lm) = Yi m (9,4>). 
The label T represents the remaining quantum numbers 
S\, S 2 , ■ ■ ■ needed to fully specify the channel. The radial 
eigenfunction R va (r) satisfies 



n d 2 + i(i + if +<s{r) 



2M dr 2 



2Mr 2 



Rva(r)=E^R va (r) 

(7) 

where, for the initial state, v labels the vibrational level. 
The spin-dipole interaction may be written as the 
scalar product of two second-rank irreducible tensors 
H sd = V p (r)T 2 ■ C 2 , where T 2 is 

ti = [si x si] I = ]T C ^ *> 2 ^< 1 - * «) S lA-» 

(8) 
C(ji,J2,j',rni,m2,m) is the Clebsch-Gordan coefficient, 
and C 2 is the second-rank tensor formed from the mod- 
ified spherical harmonics 
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The radial factor is V p (r) = b/r 3 where b = —y/6f3/h . 

To second order, the change in energy of the state |a) 
due to the perturbation H s( i is 



where 



and 



AE a = AE^+AE£K 
A£« = (a\H sd \a), 



AE^ = £ 



(a\H sd \f3)(l3H sd \a) 
p(o) p (0) 



(10) 

(11) 

(12) 



The calculation of the first order energy correction is 
straightforward. To evaluate the second order correction 
we use the method of Dalgarno and Lewis [Tfl and in- 
troduce an operator F which satisfies the inhomogeneous 
equation 



so that 

This gives 
A£i 2 ) 



[F,H ]\a)=H sd \a), (13) 



■\H Bd \a) = (E^-^){0\F\a). (14) 



(a\H sd F\a) - (a\H sd \a)(a\F\a). (15) 



As the channel states form a complete orthonormal 
basis over the angular and electronic coordinates, we can 
expand the state F\a) in terms of them: 



F|a>=X)/a'(r)|4„'>, 



(16) 



so that l|13(l becomes an inhomogeneous equation for the 
radial states f a > (r) 

(eW - h a ^j U{r) = D a , a V p (r) r" 1 R va {r) (17) 

where ($ o /|# |$ o ) = K>5 a >,a and ($ a ,\H sd \<S> a ) = 
V p {r)D a , a . 

By invoking the Wigner-Eckart theorem we find that 
the angular momentum factors in the matrix elements for 
the spin-dipole interaction between the channel states are 

D a 'a — $M s ,+m',M s +m ( _ 1) S ' S 

xC(S, 2, S'; M s , M s > - M s , M s ,) 
xC(l, 2, 1 ; to, to' — to, m) 

x(T'S'\\T 2 \\TS){l'\\C 2 \\l). (18) 

The reduced matrix elements are given by 



(t's'\\t 2 \\ts) 



d r',r osj,Si 6 s' 2 s 2 n 
x y /S 1 (S 1 + l)S 2 (S 2 + l) 



Si 


s 2 


s ' 


1 


1 


2 


5'i 


s 2 


S' 



(19) 



and 



(l'\\C 2 \\l) 
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C(l, 2, l';0, 0,0). 



(20) 



The angular momentum coefficient in Ijl9(l is related to 
the Wigner 9 — j coefficient 
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Introducing the functions g a '{r) = (r/bD a > a ) f a '(r) 
then (|17fl becomes the inhomogeneous differential equa- 
tion 
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W' + W ,y fi ( ,_ E (0) 
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III. APPLICATION TO THE v = 14 5 S+ STATE 
OF He(2 3 Si)+He(2 3 Si) 



We now consider application of the theory to the spin- 
polarized s-wave v = 14 state in the 5 E+ potential of 
two colliding metastable helium atoms. The initial un- 
perturbed state is therefore 



|a)=r- 1 i^ =14 ,a(r-)l*a), (23) 



where 



|$ a ) = |r, S = 2, M s = 2}\l = 0, m - 0). (24) 

Here T = {Si = 1, S2 = 1, f } where f denotes the quan- 
tum numbers specifying the £+ nature of the problem. 
The spin-dipole interaction coefficients l|18|) simplify to 

D a , a = 5 Ms ,+m>,2 5v i2 (-l) Ms ' 

xC(2,2,S';2,M s ,-2,M s >) 

x(r'S'||T 2 ||r2)(2||C 2 ||0). (25) 

From l|25() it is clear that only d-waves are connected 
directly to the initial s-wave state. Consequently the 
first-order energy shift for s-wave states vanishes. The re- 
duced matrix elements relevant to the second order calcu- 



lation are (r'2||T 2 ||r2) = 5r> 



r',r 



-Sr> 



r'.r 



/5/3 and (2||C 2 ||0) 




7/3, (r'0||T 2 ||r2) = 
The matrix el- 



ement (r'l||T ||T2) vanishes by symmetry of the 9 — j 



symbol. 

The energy shift for s-wave states is given to second 
order by 



f°° 1 

AE a = Y,\Da> a \ 2 b 2 R* va (r)-g a , 
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(r)dr. (26) 



Since neither R a (r) nor g a '(r) depend on the magnetic 
quantum numbers m and Ms, we may carry out the sum- 
mation to obtain finally 



AE n 
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Ka(r)^k 
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(r) 



7gr,S'=2,i=2(r)] dr. 



(27) 



IV. RESULTS AND DISCUSSION 



The differential equations JJJ for the unperturbed v = 
14 5 E+ bound state and (|17|l for the perturbed function 
g a ' (r) are solved using a discrete variable rep resentation 
and a scaled radial coordinate grid (see [ll| and Q for 
details). 

As input to the problem we require the Born- 
Oppenheimer potentials for both the 5 E+ and 1 £+ elec- 
tronic states of the metastable helium dimer. For the 5 E+ 
potential we have used the full (adiabatically and rela- 
tivistically corrected) potential Vpj (r) of Przybytek and 
Jeziorski 2] adjusted to match the experimental binding 



TABLE I: Spin-dipole induced lifetime (in units of 10 -4 s) of 
the v — 14 5 E^~ state. Results are shown for the two quintet 
molecular potentials Vp 3 (r) and Vg 5 M (r), and the two ioniza- 
tion widths TmM and TgmsM- Also shown are results for 
±20% variation in the ionization widths. Results in brackets 
correspond to a +2% variation in the singlet potential Vyi(r). 





Tm 


o.8r M 


i.2r M 


>4m 


1.578 (1.643) 
1.597 (1.630) 

Tgms 


1.779 (1.870) 
1.806 (1.850) 

0.8FgmS 


1.462 (1.510) 
1.475 (1.501) 

1.2FgmS 


V 5 


1.376 (1.429) 
1.386 (1.424) 


1.484 (1.556) 
1.499 (1.546) 


1.321 (1.359) 
1.327 (1.357) 



energy of the least bound (v = 14) state [y • In order to 
gauge the sensitivity to the potentials we have also used 
the potential Vg 5 M (r) of Starck and Meyer |13> adjusted 
to match the experimental value of the least bound state. 
For the singlet potential we use @| a potential V^ con- 
structed from the short-range Miiller et al [T3 potential 
exponentially damped onto the quintet potential at long 
range. 

In order to model Penning ionization from the sin- 
glet state we use a complex optical potential of the 
form V(r) — iT{r)/2. Two forms of the ionization width 
r(r) are used; Fm(?~) obtained from a least squares fit 
to the tabulated results in 13] and the simpler form 
IgmsM = 0.3exp(-r/1.086) advocated in 14]. The 
sensitivity of the calculation to the ionization widths was 
estimated by assigning a 20% relative uncertainty to both 
forms of the width. 

The coupling between the initial quintet state and the 
singlet state produces a complex energy shift AE a — 
AE 1 ™ — ij a /2 due to the complex form of the singlet po- 
tential. The 1/e lifetime of the least bound v = 14 state 
of the 5 E+ potential is then r = h/-f a . Using the ad- 
justed Przybytek and Jeziorski quintet potential |2| we 
obtain (see Table |TJ values of r = (1.5 78 tang) x !0" 4 



and 



(1.3761™) x 10- 



s, using the two forms 
of the ionization width 13] and 14] respectively. The 
superscripted (subscripted) numbers indicate the uncer- 
tainty due to a negative (positive) variation of 20% in the 
widths. The use of the modified Starck and Meyer quin- 
tet potential changes the lifetime (see TableHJ by < 2/zs. 
We have also studied the sensitivity of the lifetime to re- 
alistic variations in the singlet potential. We find that 
(see Table HJ a variation of +2% changes the lifetime by 
< 9 fis (a variation of —2% was rejected as it reduces 
the number of d-wave bound states from 28 to 27). We 
note that, for all cases, the real part of the energy shift 
increases the binding energy by about 10 kHz. 

Due to the disparity between our calculated value of 
around 150 /zs and the observed experimental lifetime of 
a few /xs [l[ , we considered if there were other spin-dipole 
induced processes which may contribute to the lifetime. 
It is our belief that the only possible energy-conserving 



transitions from the initial bound state are to the ion- 
ization continuum via the singlet state. For example the 
d-wave scattering states are much closer than the nearest 
lying (i-wave bound state in either the singlet or quintet 
potential, but have negligible overlap with the v = 14 
state due to the centrifugal barrier for I = 2. Due to 



the selection rule Ms< + m' = Ms + m = 2, second- 
order transitions back to s-wave states with Mg < 2 are 
forbidden. We believe that the observed lifetime is not 
spin-dipolc induced but due to some other process such as 
exchange of vibrational and translational energy through 
atom-molecule collisions. 



[1] S. Moal, M. Portier, J. Kim, J. Dugue, U. D. Rapol, 

M. Leduc, and C. Cohen- Tannoudji, Phys. Rev. Lett. [8] 

96, 023203 (2006). 
[2] M. Przybytek and B. Jeziorski, J. Chem. Phys. 123, [9] 

134315 (2005). 
[3] J. C. J. Koelemeij and M. Leduc, Eur. Phys. J. D 31, [10] 

263 (2004). 
[4] J. Leonard, M. Walhout, A. P. Mosk, T. Miiller, [11] 

M. Leduc, and C. Cohen- Tannoudji, Phys. Rev. Lett. 

91, 073203 (2003). [12] 

[5] J. Leonard, A. P. Mosk, M. Walhout, P. van der Straten, 

M. Leduc, and C. Cohen- Tannoudji, Phys. Rev. A 69, [13] 

032702 (2004). 
[6] J. Kim, U. D. Rapol, S. Moal, J. Leonard, M. Walhout, [14] 

and M. Leduc, Eur. Phys. J. D 31, 227 (2004). 
[7] P. O. Fedichev, M. W. Reynolds, U. M. Rahmanov, and 



G. V. Shlyapnikov, Phys. Rev. A 53, 1447 (1996). 

V. Venturi, I. B. Whittingham, P. J. Leo, and G. Peach, 

Phys. Rev. A 60, 4635 (1999). 

T. J. Beams, G. Peach, and I. B. Whittingham, J. Phys. 

B: At. Mol. Opt. Phys. 37, 4561 (2005). 

A. Dalgarno and J. T. Lewis, Proc. R. Soc. A 233, 70 
(1955). 

G. Peach, I. B. Whittingham, and T. J. Beams, Phys. 

Rev. A 70, 032713 (2004). 

J. Starck and W. Meyer, Chem. Phys. Lett. 225, 229 

(1994). 

M. W. Miiller, A. Merz, M.-W. Ruf, H. Hotop, W. Meyer, 

and M. Movre, Z. Physik D 21, 89 (1991) 

B. J. Garrison, W. H. Miller, and H. F. Schaefer, J. 
Chem. Phys. 59, 3193 (1973). 



